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Abstract. Up until now a complete scan in all phenomenologically relevant directions of the MSSM at the TeV scale for 
performing global fit has not been done. Given the imminent start of operation of the LHC, this is a major gap on our quest to 
discovering and understanding the physical implications of low energy supersymmetry. The main reason for this is the large 
number of parameters involved that makes it computationally extremely expensive using the traditional methods. In this talk I 
demonstrate that with advanced sampling techniques the problem is solvable. The results from the explored 24-parameter TeV 
scale MSSM (phenoMSSM) are remarkably distinct from previous studies and are independent of models for supersymmetry 
breaking and mediation mechanisms. Hence they are a more robust guide to searches for supersymmetry and dark matter. 
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PhenoMSSM 

The 105 free physical parameters of the MSSM with 
R-parity ||T[] makes a complete study of supersymmetry 
an impossible task. For practical purposes phenomenol- 
ogists had to construct models with fewer parameters at 
unification scale from which RGBs were used to obtain 
the lower-energy {Msusy) scale sparticle spectrum and 
properties. Most famous among this class of construc- 
tions is mSUGRA/CMSSM which has just 4 parameters 
and a ± sign. This has been used for providing bench- 
mark points for sparticle searches and phenomenology. 
But the approach has important limitations. Supersym- 
metry discovery definitely requires probe of large re- 
gions, in a maximal manner, of its parameter space. 
Moreover, it may be misleading if the models used to 
interpret experimental results are not realistic or not the 
most general. Hence the study of MSSM in its complete 
parameters around the electroweak scale is a more natu- 
ral approach. This will be independent of supersymme- 
try breaking models, hidden-sector physics, mediation 
mechanisms and renormaUsation group running interpo- 
lations. 

In the following sections I demonstrate this more nat- 
ural approach by applying advanced Monte Carlo tech- 
nique, called nested sampling [2] which is imple- 
mented in the currently private code MultiNest 
to fully explore the MSSM with R-parity and minimal 
flavour violation (MFV) |4] in its entire, 24-dimensional, 
most phenomenologically important parameters at the 
scale Msusy ~ 1 TeV - a set-up we call phenoMSSM |5]. 
The codes SFitter and Fittino |6] can reconstruct 
(weak-scale) MSSM parameters from collider data but 
the main goal here is to explore all the parameters, per- 
form a global fit to current indirect data and draw infer- 
ences for LHC and (future) LC physics. I start with defi- 
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nition of Bayes' theorem, describe how its variables were 
constructed for the phenoMSSM and then give sample 
results from the exercise of applying the theorem on the 
model before concluding in the last section. 

Bayesian Inference in Particle Pliysics:. Bayes' 
theorem is at the core of the algorithm used to efficiently 
explore the entire viable parameter regions of the model. 
It states that 

Pie\D,H) = P{D\e,H)P{e\H)/P{D\H) (1) 

where P{Q\H) is the parameters prior density distri- 
bution representing the conditional probability of a set 
of parameters Q given that the model or hypothesis, 
H, is true. It states what values the model parame- 
ters are expected to take. With some set of predic- 
tions (data or observables), D, obtained from the model, 
P{D\9,H) quantifies the likelihood for the model, at the 
given parameter values, to be true. For a model with 
n parameters the n-dimensional integral Z = P{D\H) = 
JP{D\9,H)P{9\H)d9 represents the evidence for the 
model. P{9\D,H) is the posterior probability density 
function and gives a measure of how well the set of pa- 
rameters predicts the given data set, D. Nested sampling 
is a general method for evaluating this n-dimensional in- 
tegral by converting it to a 1 -dimensional integral over 
a unit interval 101 ■ The sampling procedure also pro- 
duces the posterior probability distribution eqn. ([T) as 
by-product. 

PlienoMSSM parameters, 9:. The soft supersym- 
metry breaking part of the MSSM Lagrangian den- 
sity have contributions from different types of interac- 
tions, -Csoft — -^gauginos ^ -^sfennions ^ -^trilinear ^ -^higgsi 

and sources the 105 supersymmetric free parameters. 
In order to suppress CP-violation and FCNC real soft 



terms, diagonal sfermions masses and trilinear cou- 
plings, and lst/2nd generation squark masses and slep- 
ton masses degeneracies were assumed. At, Ah and 
At: are the most important trilinear couplings but we 
also include Ag = A^ because it is relevant for {g ~ 
2)fji computation f^. All the other trilinear couplings 
and neutrino masses are set to zero. This way the to- 
tal number of free parameters becomes 20. Adding to 
these the 4 most important SM "nuisance" parameters: 
{mt,mh{mh)'^^ , aeminiz)^^ , aj(mz)*^^}, makes a total of 
24 physical parameters in phenoMSSM. These are listed 
and described in Table [T] Representing them in a 24- 



TABLE 1. The 24 parameters of phenoMSSM set-up. 


Parameter 


Description 


Ml, M2, M3 


Bino, Wino and Gluino masses 




lst/2nd generation Li slepton masses 




3rd generation Li slepton mass 




lst/2nd generation sleptons masses 




3rd generation slepton mass 


ma, = mj^ = 




rricL = msi^ 


lst/2nd generation Qi squark masses 


mj^ = mj^^ 


3rd generation Qi squark masses 


= ingR 


lst/2nd generation Ur squark masses 




3rd generation Ur squark mass 


"^dR = ""'SR 


lst/2nd generation Dr squark masses 




3rd generation Dr squark mass 




top, b- and T- quark trilinear couplings 


Ae = Aji 


jj. and e trilinear couplings 




up- and down-type Higgs doublet masses 


tanjS 


scalar doublets vevs ratio 


ntt 


top quark pole mass 




b-quark mass 




electromagnetic coupling constant 




strong coupling constant 



dimensional vector 0, the combined prior for the model 
is 

n{e) = P{e\H) = 7r(0i) 7r(02) ■ • • ?r(024). (2) 

Set of Observables, D:. We use high precision 
electroweak and B-physics collider observables and the 
dark matter relic density from WMAP5 results (all 
shown in Table |2]l to study the viable parameter re- 
gions of the phenoMSSM. Its predictions for these 
observables were obtained from the 24 input param- 
eters via SOFTSUSY2.0.17 IH] for producing the 
MSSM spectrum; micrOMEGAs2 . 1 f[J\ for comput- 
ing neutralino dark matter relic density, the branching 
ratio BR{Bs ji^jj.^) and the anomalous magnetic mo- 
ment of the muon — 2)^; Superlso2.0 |14] for 
predicting the Isospin asymmetry in the decays B 
K*y and BR{b sy) with all NLO supersymmetric 
QCD and NNLO SM QCD contributions included; and 
susyPOPE ifTsIl for computing W-boson mass mw, the 
effective leptonic mixing angle variable sin^ ^eff' 



total Z-boson decay width, Tz, at two loops in the domi- 
nant MSSM parameters. These physical observables de- 

TABLE 2. A summary of the 11 observables. 
SlcoM^ error is inflated to 0.02 to accommodate the- 
oretical uncertainties. 



Observable 


Mean value 


Uncertainty 




80.398 GeV 


0.0025 GeV 




2.4952 GeV 


0.0023 GeV 


sm d^ff 


0.23149 


0.000173 


dof, X lO'O 


29.5 


8.8 


Br{h sy) X 10* 


3.55 


0.72 


nth 


1 14.4 GeV 


lower limit 


Br{B ^ ^+ li') 


5.8 X 10"** 


upper limit 




0.85 


0.11 




1.2589 


0.4758 


Ao- 


0.0375 


0.0289 




0.1143 


0.02 



rived from the model parameters form the data set, D. For 
each element. A, inDthe likelihood L(0) ^P{Di\e,H) 
was calculated. Assuming that the observables are inde- 
pendent then the combined likelihood for the model is 

L(0) (W)"''"exp[-(0,-~MOV2cT'] (3) 

i 

where (9, is the predicted value of the ith observable with 
/ = 1 , 2 , 3 , . . . , 24 and a, its corresponding standard error. 

With all the above Bayesian inference parameters 
set and ready the Multinest code was employed 
for (guided) sampling of the 24 parameters. At each 
parameter-space point the values werepassed in SUSY 
Le Houches Accord (SLHA) format fSj to the different 
particle physics software used for predicting the physical 
observables. The predictions were then checked against 
experimental values with (non)deviations quantified by 
the likelihood function. Next, the likelihoods modulate 
the parameter prior probabilities to produce the Bayesian 
evidence and posterior probability distributions for the 
model. Two prior probability density ranges, 1 TeV and 
2 TeV, were used for the purely supersymmetric param- 
eters (the first 20 listed in Table [U with the gaugino 
masses and trilinear couplings allowed to take both posi- 
tive and negative values. The slepton and squark masses 
were bounded from below at 100 GeV. The SM parame- 
ters were taken as Gaussian noise around their mean ex- 
perimental values: m, = 172.6 ±1.4, mh{mt,)^^ =4.2± 

0.07, l/a™(OTz)'*^= 127.918 ±0.018 and a,(mz)'^ = 
0.1 172 ±0.002. 

Results. Here I give some results from the explo- 
ration exercise, more results and analysis on the global 
fit to data is in progress f3\. The results are quite ro- 
bust under change of parameter prior ranges. Figure [T] 
shows that niifi is most likely around 115 to 117 GeV, 




FIGURE 1. Left: Posterior probability distribution (ppd) 
for m/jO. Right: Marginalised ppd showing chargino co- 
annihilation; colour coded: for blue to red/black for maxi- 
mum. 




FIGURE 2. ppd of gluino to neutralino mass ratio showing 
most probable phenoMSSM values around 2 and severely dis- 
favouring mSUGRA/CMSSM and AMSB. 



just above the LEP limit. The different neutralino dark 
matter (DM) annihilation mechanisms at early times of 
the universe that leads to the value of its relic density 
today are important in defining different phenomeno- 
logical regions in model parameter space. A profound 
feature of the phenoMSSM is that it shows lots of 
chargino co-annihilation as shown in Figure [1] The 
gluino to neutralino mass ratio is an interesting quan- 
tity that characterises different models of supersymme- 
try (breaking.) For instance mSUGRA(AMSB) with pre- 
dominantly bino(wino) LSP has nig/m^Q w 6(9) Olql . 
The mirage mediation ITtIi and the LARGE volume ifisll 
scenarios have the characteristic ratio less than 6 and 
between 3 to 4 respectively. Figure |2] shows the ppd of 
this quantity, providing discrimination between the phe- 
noMSSM and the other models. This feature is robust to 
prior parameter ranges. 



Conclusion:. With advanced Bayesian technique it 
is possible to fully explore weak-scale MSSM in all of its 
phenomenologically relevant parameter directions. Do- 
ing this is very important since low energy supersymme- 
try is the main focus of the LHC supersymmetry-search 
experiments. Hence such low energy approach as demon- 
strated here will be a better guide and much more realis- 
tic. The phenoMSSM has a characteristic nig /m^o < 2 



and shows the existence of chargino-neutralino (DM) 
co-annihilation. The whole procedure can be applied to 
other BSM constructions such as those with additional 
CP-violating and FCNC sources and non-zero neutrino 
masses. Moreover, using the evidence value, Z, the tech- 
nique is a powerful tool for model comparisons ifToll . 
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